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A mixed layer dynamic nodal fcr rhe structure and evolution 
cf aeicscls in a marine regime is presented. The aerosol 
spectrum is divided into continental and marine components, 
with a constant vertical continental profile assumed; both 
in and above the mixed layer. The aerosol spectrum is trans- 
formed into a reference relative humidity of 803. The 
temporal evolution of the aerosol spectrum is predicted from 
rate equations which require a specification of the surface 
producticn rate, the entrainment rate (W e ) , and the mixed 
layer depth (h) . The model was tested against the data set 
obtained in the Monterey Bay during the MAGAT 80 experiment. 
The model was initialized with both observed MAGAT data and 
an equilibrium initial value, gathered from the JASIN 
Experiment. The model was run for radii equal to 0.8, 2.0, 
5.0 and 10. C microns. The significance between the observed 
initial values and the JASIN data is that the MAGAT data 
were observed from the same air mass as the initial atmos- 
pheric data and verification aerosol data. The JASIN initial 
aerosol data are based on an equilibrium state as a function 
of only wind speed and reference relative humidity. The 
model continuously generated a correct gain or loss of 
aercscl concentrations as defined by the observed MAGAT 
data, and in most cases the model output is within one order 
cf magnitude of the observed values. 
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I. IU 1EGDUCTI0N 



There is an increased emphasis in the ability to fore- 
cast the behavior and evolution of aerosols. The scattering 
and absorption of light are influenced by aerosol distribu- 
tions in the atmospheric boundary layer and affect the 



Force is particularly interested in how aerosol extinction 
affects the use of precision guided munitions (PGM) 



PGM ' s that operate at differing wavelengths which range from 
the visible to the microwave regions. PGM's have a greater 
ability to hit a target than conventional munitions; 
however, the controlling factor is the ability of the 
guidance system to "see" the target. This ability is depen- 
dent cn the wavelength for which the sensors are designed 
and the properties in the intervening atmosphere. The 
degrading properties in the atmosphere are principally 
molecular absorption and aerosol scattering. The wavelengths 
for the guidance systems are designed such that the molec- 
ular absorption is minimised: therefore, scattering by aero- 
sols becomes the main concern once a suitable 

absorption-free window has been selected. 

The model under consideration includes the behavior of 
marine aerosol components as well as previous continental 
components. Estimating the influences of marine aerosols on 
electro-optical (EO) systems has been studied [ Barnhart and 
Streets, 1970/ . Particle sizes of interest are those asso- 

ciated with locally generated sea salts because of their 
effects on IR as well as visible wavelengths. The size 
distributions of sea salt particles show a variance of 
several orders of magnitude. 



performance of optically guided weapon systems. The Air 




The Department of Defense (DQD) has 
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The ability to forecast the behavior of aerosol extinc- 
tion from synoptic seals patterns would help in the deciding 
which type cf weapon system to employ. Because some systems 
are launched from the air, it is important that such a 
profile include the vertical distributions of aerosol carri- 
oles. Models exist for estimating vertical extinction 
profiles, but they have not been sufficiently verified. To 
do this, profiles cf acrual aerosol dara must be gathered 
and compared with the model forecasts. 

Models in current use are based on parameter iza tiens of 
the effects of relative humidity and wind speed on the equi- 
librium aerosol distributions [Wells, et al, 197?y . Recent 

evaluations have shown that these models are limited tc mean 
distr ibuticr.s (i.e., the average aerosol concentration at a 
given wind speed and humidity) /_Fairall et al, 1982a 7 . 
Models are limited because seme processes in the atmospheric 
mixed layer which affect aerosol concentrations are not 
considered, namely entrainment and inversion height changes. 

The purpose of this study is to present and evaluate a 
model which includes the meteorological processes that 
adequately describe the whole marine atmospheric boundary 
layer (MABL) . An inversion represents a cap to the vertical 
transport of surface generated aerosols, and is not 
accounted fer in previous models. The top of the boundary 
layer is capped by the marine inversion, where entrainment 
cf overlying air takes place. Because entrainment mixes 
clear (nen-marine) air into the marine layer, this process 
is as important as surface layer fluxes in determining equi- 
librium concentration. This entrainment process is included 
in the mcdel. 

Evaluation of the mcdel cutput will be dene with the data 
from an experiment entitled M ari ne Aerosol Gene rati on and 
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vicinity of Monterey Eay, Ca, during the period of 28 April 
to 9 May 1980. Ihe purpose of the experiment was to examine 
the compatibility of optical and micrometeorological propa- 
gation theory, and tc extend dynamic models of the evolving 
MABL tc include aerosol and turbulence profiles 
[fair all, 1980 and Fairall et al,198 £ ] •. Two platforms, the 
R/V ACANIA, and an aircraft were used. 

In this study, two 24-hour periods were chosen for evalu- 
ation; 3 May, beginning at 1200 PDT and 6 May, beginning at 
1800 EET. These periods were chosen because the boundary 
layer was undisturbed (no fronts closer than 100 nm) fcr 24 
hours prior to the starting times. Additionally, the times 
were selected due tc the proximity of the aircraft and 
surface ship during the experiment times. The approach was 
to describe the synoptic conditions from 24 hours in advance 
of the model forecast through the end of each forecast, and 
compare the model output with the actual MAGAT findings. 



II. background 



1. DESCEIFTIO N OF EXISTING MODELS 

Current models fcr estimating aerosol equilibrium distri- 
butions use meteorological inputs of (10-m height) wind 
speed ar.d relative humidity. These two quantities are 
considered because of their role in generation and transport 
(wind) and growth (humidity) of aerosols. The Navy's 
Wells-Munn-Katz (WKK) / Wells et al , ‘\91Tj is an example of 

this concept. The performance of this type of model has 
been studied with data obtained in the northern Atlantic and 
eastern Facific Ccean areas. The model output compared with 
these data is shewn in figure 2.1, depicting a height depen- 
dence cf tctal aerosol volume from a sample set cf eastern 
Facific data. These profiles correspond to (1) the observed 
sea salt volume, V; (2) the observed sea salt volume 
adjusted tc 80% relative humidity, V 0 ; and (3), the WMK 
predicted volume adjusted to 80% relative humidity 
(circles). it is clear from this figure that within the 
nixed layer the observed decrease of aerosol volume with 
height is less than the model predicts. The surface gener- 
ated aercscls appear to be well mixed below the inversion 
when normalized to 80% relative humidiry. 

The assertion that existing models can predict only a 
mean value appears in the results obtained in the North 
Atlantic, from the OASIN Experiment. Figure 2.2 compares a 
single radius size cf five microns and the corresponding 
model prediction. The values and trends in the predicted and 
mean results are in reasonable agreement. However, the 

standard deviation is three times the mean; if one assumes a 
normal distribution, only 67% cf the observed aerosol 
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distributions will be withir a factor of three cf -‘-he 
average. This comparison emphasizes the point than r.c manner 
how accurately a model predicts an average ae'rosci density 
at a given wind speed and relative humidity, the factcr of 
three EMS variations cannot he eliminated without consid- 
ering mere meteorological parameters. 

E. 1 EE HES BOUNDARY IAYER BODEL 

The distinguishing feature of this model is the charac- 
terization cf the MAEL, which is convectively mixed up to a 
height h, and capped by an inversion. The atmospheric 
profile is depicted in figure 2.3a, representing a cloud- 
free mixed layer where water vapor mixing ratio (g) , and 
virtual potential temperat ure (e v ) are "well-mixed", ie, 
independent of height in the mixed layer. The assumed 
vertical aerosol profile is shown in figure 2.3b. The model 
produces a 24-hcur time evclution of an aerosol spectrum, 
requiring a prediction of the following at each time step: 

1) surface production rate of marine aerosols 

2) entrainment rate at the top of the mixed layer 

3) mixed layer depth. 

Details cf how these parameters are input into the model are 
discussed in chapter IV. The mathematical relationship 
between the time rate of change of the aerosol volume spec- 
trum, dv/dr, and these three parameters is shown in the 
following equation 

dv/dt - (<wV > 0 - (w e +w km ) v)/ h 

1 2 3 (2,1) 

The model predicts the evolution of aerosol at five radii 
(0.8, 2.0, 5.0, 10.0, and 15.0 microns) of both the conti- 
nental and the marine (sea salt) components in the mixed 
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Figure 2.3 Profile of a) Mixed layer and 
b) vertical aerosol distribution. 



layer. The actual tine evolution cf each radius is cutout in 
the fern of dv/dr, with the units of um 2 /cm 3 . 

The model assumes r.c appreciable concentration of marine 
aeroscls above the mixed layer. The entrainment cf air into 
the mixed layer will not cause an element of these marine 
particles tc "escape" into the free troposphere. The ccean 
acts as a scurce for marine aerosols, primarily through the 
generation cf white caps. This is the only input for the 
marine component. The ccean also acts as a sink fer the 
continental aerosols, which are generally of a smaller 
radius. The entrainment process at the top of the mixed 
layer mixes clear (nen-marine) air into the MA3L. This 
process could be as importart as the surface layer fluxes in 
determining an aeroscl concentration. 

C. MODEL INPUTS AND AEROSOL INITIALIZATION 
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3) Frcm figure 2.3, temperature and relative humidity are 
parameterized by virtual potential temperature and water 
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vapor nixing ratio, respectively. Mixing ratios are input 
above the mixed layer, as well as the lapse rate above the 
" j u irp" . 

An initialization cf the six radii concentrations in 
three regimes is required. To do this it is necessary to 
understand the actual distributions of the continental and 
marine components in the mixed layer. For the purposes of 
evaluating the model, assumptions must be mads cn the 
distributions of the acquired MAG AT data. The aircraft 
instruments could net distinguish the chemical maks-u? of 
the individual particulates; therefore, it is not kr.cwn, for 
a given radii concentration in the mixed layer, hew much is 
continental and how much is marine. The concentrations of 
each size were calculated in the mixed layer, and then abeve 
the mixed layer. Based on the assumption that there are no 
marine concentrations above the mixed layer, we conclude 
that the concentrations calculated from above are all conti- 
nental, and therefore linerally subtracted from the concen- 
trations in the mixed layer, leaving only marine 
concentrations. 

A further adjustment of the input aerosol values is 
required, based on changes in relative humidity during the 
data collection time. We have stated that the growth of 
aeroscls is a function of relative humidity. Therefore, a 
distribution of aeroscls gathered at 90% humidity cannct be 
directly compared to another distribution gathered at 80 % 
humidity; a reference humidity is required. Consider the 
aerosol vclume spectrum 

v(r) = 4/3 Trr 3 n(r) (2,2) 
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where n (r) is the number density spectrum. V (r) 
as the volume of aerosol particles per cm 3 at 
saturation ratio of S=30%. A humidity growth 
/f airall, et al, 1982 sj is defined as 



is i s f i ns d 
a reference 
factor G (S) 



G(s) “ .81 exp(.066S/(1.58-S)) 



There are further considerations concerning advecticn. 
Clearly all aerosols do not originate locally; both above 
and below the inversion aerosols are advected into a local 
region. In terms of aerosol density, entrainment acts as an 
aerosol flux out of the boundary layer because the concen- 
trations above and below the layer are different. In the 
model, the entrainment acts on the " jump" across the inver- 
sion. For the purpose of evaluating the the evolution of a 
local concentration of aerosols, the model will neglect 
horizontal advection and further assume a negligible local 
production of the continental component. The model takes 
into consideration a Stokes gravitational fallout term, w k 
jv u, 197 s] . The fallout rates above and below the inversion 
are different because of the change in aerosol spectra, 
caused by the humidity growth factor. The Stokes velocity 
/Fairall, et al, 1962a/ is calculated from the following 
equation where p w is the density of the droplet and E is the 
kinematic viscosity of the air. 

W k = 2g(p w -p) r 2 G 2 (S) / (9E P ) 
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E. EC0I1IBFI0M AEROSOL HODEL 



As a means of comparison cf aerosol initialization 
schemes, the model offers the option of initialization with 
equilibrium values. These initial values are based on data 
collected from the JASIN Experiment £?airall, et al, 1982 J . 
A large aerosol data base was collected during equilibrium 
conditions, with respect to wind speed, during 12 hcur 
periods and then normalized to a reference relative 
humidity. Aerosol spectra were then grouped into six 
different wind speed ranges. For the equilibrium initializa- 
tion scheme in the NFS model, a pre-assigned value cf the 
five radii are used as initial aerosol values, based on the 
wind speed at the model start time. An example of the equi- 
librium aerosol spectra from the JASIN Experiment is shewn 
in figure 2.4. Note that the reference relative humidity for 
this experiment was 87 percent. The graph is referenced here 
only to shew an equilibrium type distribution; the initial 
equilibrium and MAGAT aerosol values will be shown in 
chapter IV. Once initialized with either equilibrium data or 
calculated data, the model will run in the same manner. 




Figure 2.4 Aeicsol speetpa as a function of size 
for selected wind speeds. 
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This equilibrium scheme illustrates the fundamental 
difference between the NFS model, which is a dynamic model, 
and a steady-state type model. A representative model using 
this equilibrium approach is the WMK model, mentioned in -he 
beginning of this chapter. This type of output describes the 
marine aercscl distributions as a function of wind speed 
(surface generation) , relative humidity (growth factor) , and 
elevation (vertical variation of aerosols with height 
assuming a steady state vertical transport process). The 
mathematical representation of this model /Fair all, et al, 
1982a / is 



wher e 

1 = the particle size in microns, 

u= the wind speed, 

v= 0.5 for u < 4 m/s, 

v= u -3.5 for u > 4 m/s, 

f= 1 + (v/60) 3 , 

y- 0.364-0 . 00293 v 125 , 

Z = height above sea surface, m , 
h = scale height, m (800 m for Z<1000 m) , 
a= 0.81exp (0. 06 6s/ (1 . 0 58- s) ) , 

S= H/100, (H= relative humidity in percent). 

The ether constants are given, based on the value of v. 



The equation described produces a number density spectrum, 
n(r) . This relates tc the volume spectrum V (r) , as shewn in 
equation 2. 2. 



n(r) = (r/a)* 1.6 2 (c 1 + c 2 v 5 )F*exp 




v, m/s 
v < 7 
v > 7 



Ci 

350 

0 



C 2 <5 



1000 1.15 

6900 0.29 
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III. DATA ACQUISITION AND SYNOPTIC CONDITIONS 
A. DATA ACQUISITION 

The cata were obtained from a ship and aircraft experi- 
ment, Marine Aerosol Gen era ticn and Transoort (X A G AT ) , held 
from 28 April no 9 May 1980, in a region 30 no 50 miles off 
the coast cf Monterey Bay, Ca. Data gathering was dcr.a with 
instruments mounted cn both an aircraft. and a surface ship, 
the B/V ACANIA. Aerosol data were obtained from the Airborne 
Research Associates turbo-charged Bellanca, using a Farticle 
Measuring system (PMS) Axial Symmetric Scattering Aerosol 
Probe (ASSAP) particle counter. All measured data were 
sampled every 2.5 seconds, with a two-scan average cf every 
five seconds. The scans were collected in "ladder" profiles, 
during which the aircraft made measurements at a constant 
altitude for two minutes, climbed to a new altitude, and 
repeated a new measurement run. The instrument utilized 60 
size channels from 2.8 to 14.0 micron radius. In most cases 
the ladders extended from near the sea surface (3 m) up up 
through the well-mixed boundary layer, to a few steps above 
the inversion. The elevation was generally up to 5 kilome- 
ters. A typical ladder profile contains 10 to 14 steps. The 
step heights were randomly chosen, but an attempt was made 
to keep each step height consistent between ladders. Air and 
dew-point temperatures were also measured and used tc calcu- 
late relative humidity for the correction factor mentioned 
in chapter two. Ihe aircraft also flew ascending spirals (in 
the vicinity of the ladders) during which other meteorolo- 
gical parameters were collected. This data yielded vertical 
soundings similar to those provided by radiosondes. 
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There were tw aerosol instruments fitted on the B/v 
ACANIA. One was the PMS model CSAS (classical scattering) 
and the ether was model ASAS (active scattering), controlled 
ty a IMS data acquisition system (DAS-32) with a computer 
interface. The shipboard system measured aerosols in 90 
different size channels from 0.9 to 14.0 micron radius. 

It was noted, during the flybys over the R/V ACANIA, that 
the aircraft measurements did not agree with that of the 
ship. The size distributions from the aircraft measurements 
were consistently smaller than the ship measurements, for 
radii greater than 1.C micron. In addition, the differences 
increased with radius. Since the shipboard aerosol system is 
newer and had a wider size range, and better sensitivity, 
aircraft aerosol data were corrected to agree with the ship 
aerosol data /Fair all, 1980 and Fairall et al,198_o7 . 

Profiles of virtual potential temperature and mixing 
ratio were obtained from three different sources. The reason 
for usinc these parameters instead of temperature and dew 
point, is that mixed layer inversions are more easily iden- 
tified with the former variables. The sources were the 
spiral flights from the aircraft, the radiosonde launches 
from NFS and from the R/V ACANIA. 



E. SINOETIC DATA 



urface and 


500 
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Hite image 
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synoptic 
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weather maps. In addition to data collected from the R/V 
ACANIA and the aircraft, local weather data were also avail- 
able from the U. S. Army's Fritzsche Field weather facility. 
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Figure 3.1 



Surface and 500 mb anaylsis for the 
at 0500 PDT 2 , 3 and 4 Hay 1980. 



western 0. 
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Figure 3.2 



GOES West satellite imagery. 



0915 PDT 2 Bay, 



1980. 
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Figure 3.3 Same as Figure 3.2 except 1245 PDT , 3 May. 
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Figure 3.4 Same as Figure 3.2 except 1645 PDT , 4 May. 
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Figure 3.5 Same as Figure 3.1 except 5, 6, and 7 Bay, 1980 
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Figure 3.6 



Same as figure 



3,2 except 1745 PDI , 



5 May. 
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Figure 3.7 



Same as figure 



3.2 except 1645 PDT , 



6 May. 
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Figure 3.8 Same as Figure 3.2 except 1745 
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Fort Crd, Ca. Satellite photographs and NOAA maps for the 
two different time periods (2-4 May and 6-7 May) are shewn 
in figures 3.1 through 3.8. Each three- day sequence depicts 
the 24 hour period being considered; 24 hours before and 
after. Atmospheric soundings for the two model runs are 
shown in chapter IV, figures 4.1 and 4.2. 

C. SINOETIC CONDITIONS 

Several weak frontal systems passed through the area 
during the experiment. Showers occurred during the first and 
last cays, associated with the fronts. Low cloudiness and 
fog cccured during the morning from 29 April to 5 May, with 
fog returning again cn the last day. 

At the beginning cf the period, the area was dominated by 
a slowly eastward migration of a cut-off low at 500 mb. By 
early morning on 2 May, the area was under the influence of 
a weak ridge. On 3 May the area was under divergent flew at 
the upper level. An upper level low had formed off the Baja 
Ca. coast cn 4 May, leaving the area under the influence of 
a col. Cn 5 May the area was between an upper-level trough 
and ridge, and by 6 Kay the area was on the back side of the 
trough. Eecause cf the deepening of the trough, the area was 
still cn the back side of the trough on 7 May. A new upper 
level trough formed end approached the area on the final two 
days cf the experiment. 

Surface winds were relatively light, 0 to 10 kts, at the 
beginning cf the peried, and increased toward the end of the 
peried tc 16 kts, gusting to 22 kts. 

An important feature of these interpretations is the 
nature of the mixed-layer, eften topped by an inversion, 
with reguard to stability, and therefore mixing intensities. 
It is assumed that tie mixing becomes greater as conditions 
become mere unstable. 
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At the start cf the experiment, the mixed layer exhibited 
stable tc slightly unstable conditions, until about 1800 
local time (PDT) 28 April, when conditions became rncr^ 
neutral. The neutral condition remains until 1 May when 
conditions once again become stable. A weak frontal passage 
before 0500 PDT on 29 April dees not appear tc affect the 
mixed layer profile. The layer remains stable until a 
frontal passage on 2 May when conditions become neutral and 
remains sc until 5 May. On the morning of 5 May, conditions 
are slightly stable, but returned to neutral on 6 May, 
despite a frontal passage at 1300 PDT on 5 May, and 



remaining neutral to the end o 
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1200 TIME (LOCAL) 1200 

05/03 05/04 

Figure 3.9 Plots of wind ?pe§d, air and sea surface 
temperature, and BH beginning 1200 PDT 3 May. 
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RELATIVE HUMIDITY OO TEMPERATURE (C) WIND SPEED (M/S) 






figure 3.10 Same as figure 3.9, except beginning 1800, 6 May. 
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Figure 3.12 Same as figure 3,11 except afternoon of 3 Hay, 1980. 
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Figure 3.16 Same as figure 3.13 except afternoon of 7 May 1980. 




Figure 3.17 Cruise track of R/v 

42 



ACANIA 3 May 1980 



123 ° 




Figure 3.18 



Saae as figure 3.16 except 4 May 1980. 
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Figure 3.19 
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Same as figure 3.16 except 6 and 7 Hay 1980. 
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were chosen; 3-4 Hay and 6-7 Hay. Figures 3.11 through 3.16 
indicate the routing of the aircraft flights for the two 
time blocks, and figures 3. 17 through 3.19 depict the ccurse 
cf the E/V ACA HI A during the same times. 

Atmospheric data from the aircraft soundings were used 
instead cf the radicscndes frcm the R/V ACANIA. This was 
done fcr twc reasons: first because the aircraft soundings 
contained mere levels, ar.d second they were in closer prox- 
imity tc the aircraft aeroscl ladders. Atmospheric and 
aeroscl cata for the first run was taken from flight 6, on 
the morning of 3 May. The locations and times are shewn in 
figure 3.9. The aircraft was approximately 45 miles frcm the 
coast when a sounding was taken at 1143 PDT, and an aerosol 
ladder (111) was taken at 1200 PDT. Data for the two verifi- 
cations of atmospheric parameters and aerosol values were 
chosen based on their proximity to the initial data. A 
sounding and aeroscl ladder on the afternoon of 3 May 
(sounding time of 17 14 PDT, ar.d aerosol ladder 113) approxi- 
mately eight hours frcm the initial time was chosen for the 
first verif icaticn, and is depicted in the route cf flight 
1, figure 3.10. Data for the second verification was chosen 
from flight 8, figure 3.11; sounding time 1151 PDT, and 
aeroscl ladder 116. Cata for the two verifications are beth 
within 2C miles cf the initial data. 

Data for the second model run is shown in figures 3-12 
through 3-14. The model was initialized with data frcm the 
1-732 PDT sounding and ladder 122, approximately 25 miles eff 
the Mcnterey coast. The first verification data was taken 
from flight 12 the next morning approximately 17 hours 
later; sounding time 1743 PDT, and ladder 123. The seccnd 
verification data was taken from flighx 13, approximately 24 
hours frcm the initial data; sounding time 1836 PDT and 
ladder 124. Data collected for the two verifications were 
within 5 miles of the location of the initial data. 
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Additional wind speed, direction, ar.d 3ST collected from 
the E/V ACANIA were used for initialization and verification 
cf the two model runs. As mentioned in the model descrip- 
tion, up to 10 forecast values for each can be input into 
the medal. However, due to the continuous movement of the 
E/V ACANIA during the MAGAT experiment, where SST values 
were charging because of strong coastal gradients, a single 
value cf SSI was used. This value was chosen when the S/V 
ACANIA was closest tc the areas where the aircraft data was 
taken fcr each cf the two initialization sites. The wind 
speed and direction cata frcm the E/V ACANIA was considered 
representative of tlie data collection sites, and therefore 
input into the mcdel runs at approximately 3 hour intervals. 



> 
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17. MODEL INITIALIZATION AND RESULTS 

A. DIGITIZED ATMOSPHERIC SOUNDINGS 

The atmospheric irputs for the model (chapter 2) are: 

1. The mixed layer equivalent potential temperature (C) 
anc specific humidity (gm/kg). 

2. The "jump" discontinuities in each of the above values. 
The jump meaning the difference between the mixed layer 
value and the value at the top of the inversion. 

3. The lapse rate for each parameter value atcve the 
inversion . 

4. The depth of the well mixed-layer. 

To simplify the process of calculating these values from the 
soundings, a digitizing scheme was designed by the 
Environmental Physics Group at the Naval Postgraduate 
Schocl. This scheme transforms the sounding into the struc- 
ture that describes the necessary model inputs. Examples are 
shown in figures 4.1 and 4.2. Note that the "jumps" are 
depicted as occuring in an infinitely thin layer. 

E. ADDITIONAL MODEL INPUT AND ADJUSTMENTS TO THE DIGITIZED 
SOUNDINGS 

The reason for making adjustments to the atmospheric part 
cf the model is to provide a better basis for evaluating the 
behavior of the aerosol prediction. This is optimized when 
the model is producing the best possible forecast cf the 
atmospheric parameters. The parameters are forecast very 
well in both model runs, with these adjustments mads. The 
model is accurately describing the height of the inversion 
through adjustments in the subsidence rate, and the lifting 
condensation level (LCL) is accurately generating clouds 
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through adjustments tc the temperature and specific humidity 
which are within the limits of instrumentation error. The 
cloud verifications were made with satellite data and obser- 
vations from the R/V SCANIA leg. 

An initial input for the model is the subsidence rate. 
Using a first guess cf -.0 05 m/sec, the model generated a 
plot cf the inversion height for the 24-hour perrod. Eased 
on the verification data, the subsidence rate was further 
adjusted tc bring both the predicted and observed inversion 
height values together. Then a small adjustment was made to 
the mixed layer values to match the stability and cloud/ 
cloud free patterns with what was actually observed during 
the period. In the first model run (3 May) , the potential 
temperature was increased 0.5 degrees and ~he specific 
humidity was decreased 1 gm/kg. The resultant subsidence 
rate necessary tc match the predicted and observed inversion 
height was -.0042 m/sec. These adjustments yield an atmos- 
pheric profile similar to what was actually observed. The 
air mass became slowly saturated enough to generate clouds 
in the upper part of the mixed layer at about 0300 PD T on 4 
May. Flcts of the inversion and LCL behavior, and the temp- 
erature and humidity profiles are shown in figure 4.3. In 

i 

the plot for LCL and inversion height, the letter "H" 
depicts the verification times and values for the inversion 
height, and the letter "L" similarly indicates values for 
the LCL verifications. The letters "T M and "Q" in the temp- 
erature and specific humidity graphs indicate verification 
values fer each respective value. 

The atmcspheric plots for the second model run (6 May) 
are seen in figure 4.4. For this run the specific humidity 
was reduced from 6.48 to 6.3 gm/kg, the potential tempera- 
ture was unchanged, and a subsidence rate of -0.01 m/sec was 
input for the run. 
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Figure 4.1 



Digitized sounding for first model run. 
3 May 1980, 1143 PDT. 




Figure 4.2 



Same as figure 4. 2 . except for second run, 
6 Hay 1980, 1749 PDT. 
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Figure 4.3 Plots of inversion and LCL heights, temperature, 
and humidity profiles, first run. 
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Figare 4.4 Same as Figure 4.3, except for second run. 
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C. AEROSOL IHIT IALI2 ATION EEOCECORES 

The aerosol data were gathered in ladders, usually up to 
five kilometers, and were of actual concentrations of parti- 
cles for a given increment cf radius, n (r) . This number 
contained bcth marine and continental components, collected 
at the ambient relative humidity. 

The initialize ticn scheme was performed in three main 
steps. The first step was to identify the marine and conti- 
nental components. By using the digitized soundings for each 
model run, a determination was made of two height values: 
cne representing the well mixed-layer and the other repre- 
senting the free troposphere, above the inversion. Eoth the 
mixed layer and free troposphere aerosol counts for each of 
the five radii were calculated. The second step was to apply 
the humidity growth factor to the values and calculate 
dv/dr . The relationship between n (r) , which was measured, 
and dv/dr is given by equation 2-2. The final step was to 
subtract, fcr each racius, the free troposphere dv/dr value 
from the mixed layer dv/dr value. With the assumption of a 
constant prefile of continental aerosols through both the 
mixed layer and into the free troposphere, this subtraction 
leaves in only the mixed layer marine components, normalized 
to 80 percent relative humidity. The values of dv/dr used 
to initialize the mcdel, both equilibrium and MAGAT, are 
shown in tables 4.1 and 4.2. 

This three step prccess was dene for both model runs, and 
for twe verification times within each run. As described in 
chapter III, the verification ladders were chosen to be as 
close as possible tc the position of the initial aerosol 
data. In both cases, the verification data were within 15 
miles cf the initial data. The plots of the aerosol output, 



53 



tcth for the MAGAT and equilibrium initializations, is 
depicted as loglO (dv/dr) against the 24 hour time period, 
in half hour time steps. In both (initialization) cases, the 
"X" values represent the MAGAT verification values. In both 
model runs, aeroscl output for the MAGAT initialization 
scheme is represented by a continuous line, and the dashed 
line depicts the run for the equilibrium initialized cases. 
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Whet considering the equilibrium initial values from 
table one, it is easy to see the relationship between wind 
speed and production. In almost all cases, the aerosol 
volume counts increase as the wind speed increases. The 
distribution of radii for a given wind speed, however, is 
not as well defined. For lower wind speeds, the larger 
radii are few ir nunber, and as the wind speed increases 
there is a dramatic increase in the number of larger parti- 
cles. This is reasonable due to white cap production. Also, 
there appears to be a slight maximum of particles at 2.0 
microns, reguardless of wind speed. The 15 micron initiali- 
zation values in the tables and the model output for the 
same si 2 e are not shown. This is because the observed 
initial values collected from the aircraft were zero. 

In comparing the MAGAT initial values, one finds that the 
smaller radii are approximately twice the equilibrium count 
for even the largest wind speed. It is interesting to net® 
that for such a large variation in initial wind speed for 
the two cays, there is very little variation in the initial 
values, for all the the particle sizes. Based on results of 
Monahan /et al 1983/ , it was also surprising to find no 
volume counts for the larger particle sizes since wind speed 
was above 10 m/sec. For the strong wind case, there were no 
aerosols above 10 microns collected from the aircraft. 

E. J-IBSl MODEL BOH EESOLTS 

The plots of dv/dr over the 24 hour period for the first 
model run are shewn in figures 4.5 and 4.6. For each radius, 
the solid line indicates the profile with the MAGAT initial- 
ized values, and the dashed line represents the profile with 
the equilibrium initialization. 
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For the three s mailer sizes, the model sho'«s little 
change ever the period, however the model does produce a 
trend in aerosol production. Whatever volumes were used to 
initialize, the number tends to stay within one erder of 
magnitude ever the period. The verification values compare 
very well to the model run with the MAGAT initialization. 
The model and verification dv/dr values differ by less than 
10 units. Eor the larger size, 10 microns, note the guick 
production of aerosol when initialized with with a zero 
value. Within four tours, the model is producing close to 
the same values as the equilibrium case. From this feur-heur 
point, the two cases behave very much the same. There is a 
diurnal variation evident in the larger size cases; a lag is 
noticeable in the production through the night time, with as 
increase after sunrise and through the morning. The verifi- 
cation values for these larger sizes are not good for either 
verification scheme. 
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Figure 4.5 Aerosol plots for first run, with a) r= 0.8 ua, 

and b) r= 2.0 um. 
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Figure 4.6 Same as figure 4.5, except a) r= 5.0 um, 

and 6) r= 10 um. 
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E. SECOND MODEL RUN EESULTS 
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Figure 4.7 Aerosol plots for second run, a) r= 0.8 um, 

b) r= 2.0 um.. 
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LOG 1 0 dV/dR, ( um A 2/cm A 3 ) LOG10 dV/dR, (um^/cm^) 
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For RRDIUS= 10 um 
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Figure 4.8 Same as figure 4.7, except a) r= 5.0 um and 

b) r= 10 um. 
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S0MM2RI OF RESULTS 



The mcdel is producing, with a minimum of adjustment, an 
accurats description of the atmospheric boundary layer 
parameters. Once a subsidence rate is selected such that the 
24-hcur predicted inversion height approaches the observed 
value, the values cf temperature and humidity only require a 
minimum cf adjustment. In both cases the verification values 
cf temperature ware within 2 * C, and the values of specific 
humidity within 1 gm/kg. More significantly, the trends 
ever the 24-hour periods fer both temperature and humidity 
were accurately described by the model, as defined by the 
verification data. 

With respect to the aeroscl input, the model is also 
producing the correct trend in behavior. There is evidence 
in the output fer the trend in aerosol production to be a 
function of not cnly wind speed and relative humidity, but 
also the concentre ticn of the initial aerosol values. This 
is seen in the first run, particle sizes 0.8 and 2.0 
micrcrs. The atmospheric parameters are the same in beth 
cases, fer each radius; but it is the initialization scheme 
that generates a different trend in each case. The mcdel run 
with the MAGAT initial values produces a loss in production, 
and the verification values are within one order of magni- 
tude cf the model values. With the equilibrium case, the 
trend is an increase in production over the time period, 
with the mcdel verification values two orders of magnitude 
away from the observed verification values. 

The production trends in the second model run, with both 
initialization schemes, are much the same. In this case 
there is less of a difference between the model verification 
values fer the two initializations. With the exception of 
the smallest radius, 0.8 microns, the verification values 



for all the radii of the second run are within one order of 
magnitude of the model output. The fact that the MAGAT 
initial values are closer to an equilibrium value at a 
higher wind speed seems to imply a less sensitive model, 
when initialized with a higher wind speed. Accordingly, 
aerosol production is high;- and the model accurately 
predicts the behavior of the aerosol in the stronger wind 
case. 
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